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THE DYNAMICAL INSTABILITY OF STATIC, SPHERICALLY 
SYMMETRIC SOLUTIONS IN NONSYMMETRIC GRAVITATIONAL 

THEORIES 

M. A. CLAYTON, L. DEMOPOULOS, AND J. LEGARE 

Abstract. We consider the dynamical stability of a class of static, spherically-symmetric 
solutions of the nonsymmetric gravitational theory. We numerically reproduce the Wyman 
solution and generate new solutions for the case where the theory has a nontrivial fundamental 
length scale ■ By considering spherically symmetric perturbations of these solutions we 
show that the Wyman solutions are generically unstable. 



1. Introduction 

The goal of this work is to understand the role of a class of static, sphericany-symmetric 
solutions of the nonsymmetric gravitational theory p|, ||] (NGT), which are also solutions to 
Einstein's unified field theory Q (UFT). In Section |2| we show that the time-dependent field 
equations for this sector are equivalent to an E-KG system with metric g and a scalar field 
■0 with a positive semi-definite (coordinate-dependent) self-interaction potential. This result 
allows us to use a scaling argument to show in Section ^ that there are neither solutions with 
nontrivial ip and globally regular g, nor solutions with an event horizon and nonzero tp in 
the exterior (no hair). The remaining class of static solutions (regular and no event horizons 
outside of the origin) consist of the Wyman solution ^ and generalizations to solutions 
of the theory that possesses a non-zero inverse length scale /i, both of which are numerically 



generated in Section 3.1 



The issue at hand is whether these solutions may be considered as a possible endpoint of 
gravitational collapse in NGT, a matter which has been discussed in the literature |^, |7|, ^, 
p. To answer this question we investigate the dynamical stability of these solutions under 



spherically symmetric perturbations. Employing both a variational argument in Section 4.1 
as well as determining the ground state and its eigenvalue directly, in Section ^ we find that 
these solutions are dynamically unstable. The instability of the = case is analogous to that 
of an Einstein-Klein Gordon system (E-KG) where it was found that the static solutions with 
nontrivial scalar field (also attributed to Wyman |l^) are unstable [11|. 

2. The Wyman Sector Field Equations 
Here we will give the (coordinate frame) action for NGT as presented in |12] 

S^cT = J A V=g(-g'^"^M^ - s'-^'^l^.iWU + lo^E^^'W^W, + FA, + i/i^gl^'^lgiH) • (2.1) 

The inverse of the fundamental tensor is defined by g^'^g^^ = gwyg^^ = Su, and g := det[g^,y]. 
The covariant derivative is characterized by the torsion-free O^f,^] = 0) connection coefficients 
V^[3jy] = Tjiud'y, not required to be compatible with any tensor. What are commonly treated 
as the antisymmetric components of the connection coefficients are considered as an additional 
tensor Aj^u. 
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The Ricci-like tensor that appears in the action is split into two contributions = R^y + 
R^y] the first is identified as the Ricci tensor defined from Vj^y (and reduces to the GR Ricci 
tensor in the limit of vanishing antisymmetric sector), and the second contains contributions 
from the antisymmetric tensor field KJiu (A^ := A/^-y). These are: 

R,. =d,[TZ,] - Id^m^] - \d,[TZ^] + Ti^T], - r^,rj^, (2.2a) 
< =v^[a]t;, + V[^[A],] + a;^,a^,^. (2.2b) 

In addition to the vanishing of the torsion tensor, the compatibility conditions of NGT are 



A^ = 0, V,[^g][H = aV^g(^^)H^„ = iag(^^)VF„ (2.3a) 

and the field equations are given by 

= R^l + - hW^Wy - Ifi^M^y = 0, (2.4) 

where 

= g[H - Mi^sg^^^^ + Ig.^gt^'^lgM. (2.5) 

The action ( |2.1[ ) and field equations (|2.3[) and ( p.4| ) encompass those of the 'massive' theory 
P, HI, when a = 3/4, 'old' NGTg with vanishing source current for a = and ^ = 
(which is equivalent to UFT |3|), and recovers GR in the limit that all antisymmetric 
components of the fundamental tensor are set to zero [12, 

The general form of the spherically-symmetric fundamental tensor [ p!^ consists of the general 
form of a spherically symmetric metric with the additional antisymmetric components g[oi] = 
oo{t,r) and g[23] = f{t,r)sm{0); in this work we consider the Wyman sector |Q, |5| defined 
by choosing gjgi] = 0. The remaining components will be parameterised by the symmetric 
sector functions and A, and a dimensionless field ijj parameterising the mixing of the angular 
components 

\Sifiu)\ = diag(e'',-e^,-r2cos(V'),-r^cos(^)sin2(6')), g[23] = sin(V') sin(6'). (2.6) 

From this form of the fundamental tensor one finds the following combinations of the field 
equations (|2.4| ) which will prove useful: 

noi=}{dt[X]-^dtmdrm}, (2.7a) 

e-7^oo + e-^7^ll =ie-^{a,[i/ + A] - Le^-{dtmf " fl^r-lV'])'}, (2.7b) 

^(cos(V)7^22 + sin(V)7^23/sin(0)) = 



2r 



1 e-^{dr[i^ - A] + f (1 - e^ cos(V')) + |^Ve^sin2(V')}, (2.7c) 



-;^(sin(V')7^22 - cos(V)7^23/sin(0)) =le-^{-e^-d![i,] + le^-''dt[iy - X^] 

+ 8^14^] + ^dr[iy - X]^r[^P] + Idri^P] - ^e^sm{^|;) - ieVsin(2V)}. (2.7d) 

Examination of the field equations (|2.7l ) reveals a strong similarity to an E-KG system. 
Pursuing this analogy we introduce the Riemannian metric 

g = diag(e^, -e^, -r^, -r^ sm^{e)). (2.8) 

Denoting the Ricci tensor derived from g as R^u, the field equations may be written in the 
E-KG form: 

Roo =U9tm)^ - Ifi^e'sm^iP), Roi = '^dt^dM, 

Ru =U9r[^]f + i/^'e^ sin2(V;), i222 = 1 - cos(V) + sin2(V;). 
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Additionally we have the wave equation ( |2.7d ) for -0, which may be written as 

r^v^v.iV'i+^^^m = 0, 



(2.10) 

(2.11) 

to date 

the analogy with an Einstein-Klein Gordon system has not been realised. This identification is 
useful since the E-KG system is known to possess a unique light cone determined by the metric 
g, and allows us to identify the metric g (up to a confirmal factor) as a physical measure of 
spacetime. This is a nontrivial statement in the UFT in general given the existence of multiple 
It is also straightforward to show that if ([ 



where <5^[V'] is the functional derivative with respect to ip of the potential defined by 

V:= ^(l-cos(V)) +i/x2sin2(0). 
Although the field equations for this sector have appeared in various places In 



20 



light cones |21 



are written as G 



It 



then the Bianchi identity Vjy[G]'^^ = is satisfied by T^,^ by virtue of the wave equation ( |2.10| ). 
Thus although the Bianchi identities for the nonsymmetric theory |13| do not in general imply 
separate conservation laws for a general relativistic sector and an antisymmetric sector written 
as a simple matter field, for the Wyman sector this is realised. 

This provides some motivation for studying this sector; considered on its own it suffers from 
neither the bad asymptotic behaviour noted in |22, ^] nor the linearization instability discussed 
in and is reproduced by many of the nonsymmetric field actions. Thus one conjectures 
that when an acceptable nonsymmetric action is found, it will also reproduce this sector. We 
show that if we identify g with the metric of an E-KG system, the static, spherically-symmetric 
Wyman spacetimes with nonvanishing tp and the E-KG spacetimes with nonvanishing scalar 
field (also attributed to Wyman |1C, [23]) have similar properties. Thus we expect that they 



have similar dynamical stability properties [T^, 26|, a result that we will establish below. 



3. The Static Background Spacetimes 

We now turn to the task of generating the static background fields i'q, Aq, and -00 • Since we 
are going to show that these solutions are dynamically unstable, it is worthwhile to include a 
proof that the solutions for which g is globally regular or has an event horizon with regular 
exterior consist precisely of the Minkowski spacetime and the Schwarzschild solution. The 
argument is an extension of that given by Heusler |2^ . 

Writing e'^ = n^(r)(l — 2m(r)/r) and e"** = 1/(1 — 2m(r)/r) and inserting the fundamental 
tensor ( |2.6| ) into the Lagrangian density ( |2.1| ) we find the effective action for static configura- 
tions (using J d^l = 4n) 

4a, H - -r^(l - —) {drii'])'^ - 2(1 - cos(V)) - -r^^ sm\i;) 



eff 



dr n 



(3.1) 



2 V r ' 2 

It is straightforward to show that variations of this action with respect to n(r), m{r) and ^{r 
are equivalent to the static limit of (p.Tj). 

Assuming that we have an asymptotically flat solution to the static field equations 

m ''~'°°> Ms + o(l/r), n ''~'°°> l -|- o(l/r), ip ''^°°> o{l/r), 

with either a regular centre at r = 6 = so that 

mor -|- ©(r^) with uiq < 1/2 



lim ih is finite, 

r-»0 



lim n is finite. 



m ' 



(3.2) 
(3.3a) 



or an event horizon at r = 6 where 



lim ip is finite, lim n is finite, 

r— »b r— >fe 



lim m = 6/2, 



(3.3b) 



then the one-parameter family of fields nx, mx, and ipx^ where for example mx{r) 



m 



((1 



A)6-|- Ar), have the same boundary values. The action S'^ := S^^[nx,mx,ipx] must therefore 
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have a critical point at A = 1, and so c?;^ [S"^^] |a=i = 0. Performing this variation on ( |3.1| ) 
results in 



oo 



v2 



dx[Sf] 1^^^ = 47r / dm i(r(r - b) - b{r - 2m)) [drlip])' 

J b 

+ 2(l-cos(V')) +|(r^ + 2r(r-6))/i2sin2(V') =0. (3.4) 

The final two terms are manifestly positive semi-definite, and the first may be seen to be as 
well by the following argument: the coefficient of (9.r[V'])^ vanishes at r = 6 for either a regular 
centre or a horizon, and its derivative 2(r — 6) + 2bdr[m] is non- negative for r > b (this follows 
from the variation of ( |3.1| ) with respect to n(r), which shows that > 0). Therefore all 

terms in (^^) must vanish separately and we have shown that no solution exists other than 
the trivial one '0 = 0. That these results are a straightforward extension of those for an E-KG 
system as derived by Heusler is due to the positive semi-definiteness of the potential ( 2.11| ). 



We therefore know that static solutions to the field equations that have nontrivial antisym- 
metric components of the fundamental tensor must have neither event horizon nor be globally 
regular. It is to the numerical determination of such solutions that we turn next. 

3.1. Numerical Results. Throughout the remainder of this work we will make use of the 
dimensionless radial coordinate scaled by the asymptotic Schwarzschild mass parameter Mg of 
the system: x := r/{2M^). In the case of a non-zero length scale we will also make use of 
this coordinate, additionally defining Ji := 2Ms^. We then re-write the static limit of ( ^.7| ) in 
a form appropriate for the numerical integration implemented below: 

9x[e^°] =e^°{f(5,[V'o])' + ^(l-e^«cos(0o)) + i/i'xe^°sin2(0o)}, (3.5b) 

dliM + Ml + cos(Vo) - i/i'x'e^" sm\^Po))^,[^Po] 

- ^e^o sin(V'o) - ^/i'e^« sin(20o) = 0. (3.5c) 

We will employ the simple shooting method |^, Section 7.3.1] on ipQ to numerically generate 
a solution to (|3.5|). Initial values for e^^^, e^°, ipQ and ^^^[^o] a-re given at a small x and ( p.5| ) are 



integrated outward using a fourth-order Runge-Kutta [28, Section 7.2.1] routine developed by 
the authors. It is therefore necessary to have the small-x behaviour of these solutions, which 
may be determined from ( |3.5| ). Assuming that both of e'^" and e"**" vanish as some positive 
power of x as X — > 0, we find 



Fx'J^^^~'\ e^o ^ ^o~-^ln(x.)+5, (3.6) 

(cp-l)^ ap-1 



independent of the value of jl. In the /i = case we know from ( A. 101 ) that F = G = 1 and 



from dAl^ ) that B = tan-^(-s). As we show in Figure |l|, the numerical solutions reproduce 



the analytic results quite well. By shooting on B we have checked that the shooting algorithm 
reproduces B = tan~^(— s) to machine accuracy. 

The large-x behaviour of the solutions is determined by requiring that the metric components 
behave asymptotically like the Schwarzschild solution: 

e'^o-l-i e^»~l + i. (3.7) 



Inserting this into (3.5c) we find (keeping only the asymptotically dominant terms and noting 
that we require that -0—^0) 

9l[M + ld,[M - 4Vo - /iVo ~ 0. (3.8) 
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From this the asymptotic form of ipo is determined: 



A/x'^ for /i = 
Ae'i^'^/x for /i / 



(3.9) 



and from ( A.8D we know that for /i = we must find A = s/12, a result that has been verified 
numerically. 
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Figure 1. The Wyman solution (/i = 0): the solid line corresponds to s = 0.2, 
the long-dashed line to s = 1, the dotted line to s = 10, and the short-dashed 
line to s = 100. The first six analytic points from (A. 13) are indicated by 
diamonds on the plot of "00 • 



For the /i 7^ solutions things become more complicated. Although the constant F may 
be freely chosen by rescaling the time variable, we still must determine the two constants G 
and B, which may be shot for by requiring both -00 — > and (|3.7| ). In order to avoid shooting 
on two variables we proceed as follows: We choose a value of s and /i and assign F = G = 1 
initially, then shoot on ipQ to determine B. This will result in a solution with asymptotic 
behaviour e^" ~ (1 -|- a/x) and e''" ~ /o(l — a/x). Re-scaling a; by a; ^ ax and renormalizing 
e'^" results in a solution with the correct asymptotic form for the metric functions, while the 
a; ~ forms become 



/o-^(ax,)2/('^p^i) e^° ~ 4Vl + s^{aXsf^''/^^''~^'^/{ap - if 



(3.10) 
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The constant B determined by shooting is shifted hy B ^ B — 2am, ln(a)/(ap — 1) and, most 
importantly, we have fl — > aji. 

In Figure we show the numericahy generated solutions for various values of s and fi. 
We find that initially for increasing values of /2 input into the above rescaling procedure the 
re-scaled /i is also increasing, however it then turns over and decreases towards zero. 




012345 02468 10 



X X 



(a) s = 10: the solid line corresponds to ji ~ 
0.01001, the dashed line to /x « 1.551, and the 
dotted line to the ji ~ 2.418 solution. 



(b) s — 100: the solid line corresponds to /i ~ 
0.009995, the dashed line to /i « 0.1471, and the 
dotted hne to the jl ~ 0.1625 solution. 



Figure 2. A selection of solutions with /i / 0. 



4. Dynamical Stability 

Finally we turn to the investigation of the dynamical stability of these spacetimes, i.e., 
whether they are unstable against linear perturbations. Two different methods will be em- 
ployed: a variational approach which casts the equation for the perturbative modes as the 
solution of a stationary Schrodinger equation and uses a one-parameter family of wave func- 
tions to put an upper bound on the lowest eigenvalue, and a shooting method to numerically 
generate both the lowest energy mode as well as the eigenvalue. 

In either case we make use of the single Fourier- mode expansions of the perturbing fields ui, 
Ai and ipi: 

^{tjx) = vq^x) + uii^x) cos{LJt), A(t, x) = Ao(x) + Ai(x) cos(a;t), (4-1) 
ip{t,x) = ipoix) + lipi{x)cos{u)t). 
The first-order perturbation equations from ( p. 7] ) for the symmetric functions give 

Ai = ^dM^i, dM = -h^llMi'i - IdM^i + Id^mdM (4.2a) 
and the equation for ipi 

dlli^i] + id,[uo - Aol^.iV'i] - e^^-^^oyfxlVi + cDVo-'^>i = 0, (4.2b) 

where 

V[x] := e'^^-^^I^^^.K - Ao] - U^Mf " j{dMfd,[,^o - Ao] 

+ f e^o smii^o)dM + :^e^« cos(Vo) + /i'e^° cos(2V'o) + i/i^xe^" sin(2V^o)5:.[^o]}. (4.2c) 
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From the small- a; behaviour of the background fields (p.6p we find 



(4.3) 



from which we have that ipi ~ x{A + B\n{x)). We will only consider cases where B = 
(corresponding to r = below) and the perturbation is therefore finite at x = 0. Once the 



(arbitrary) amplitude A is chosen, ( [1.2b ) is integrated numerically, varying the value of (Z"^ and 



requiring that V'l ^ as x ^ oo. In Figure ^ we show the results of this procedure, giving 
the value of o)^ for the jl = solutions, as well as a selection of /i / solutions. In addition, 
we were unable to find any additional unstable modes. 




-0.05 



-0.15 




(a) jl = 0: The solid line represents the en- 
ergy of the unstable mode determine d nu meri- 
cally from the perturbation equation (4.2), and 
the shaded region the upper bound determined 
from the variational method (If. 20 



(b) jl ^ 0: The energy of the unstable mode for 
(bottom to top) s = 20, 23, 25, 30, 40, 50, 70, 
100 



Figure 3. Unstable mode energy for ^ = and ^ 7^ solutions respectively. 



4.1. Variational Results. Following Jetzer and Scialom [11| we introduce a radial coordinate 
p that eliminates the second term in ( [4. 2b ) 



P 



dxe"2(^o-^o), (4.4a) 



(4.4b) 



chosen so that p = at x = 0. From ( p.6D and (|3.7D respectively we find that 

[ JG/Fx'^/{ap - 1) as X ^ 
I X — ln(x) + constant as x — > 00 ' 

p is determined numerically using the x — > result and integrating outward. 
This coordinate transforms ( 4.2b| ) to the stationary Schrodinger-like equation 

-^ + ^[x(p)]^i=cI^Vi. (4.5) 

We will consider perturbations in the Hilbert space C^[dp, (0, 00)), i.e., with the norm 

r°° 1 

dp(l){p)(l)'{p) = (ixe-2('^o-^«)</.[y9(x)]0'[p(x)]. (4.6) 
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From the asymptotic and small x forms given earlier, we find 



2//)2 for /i = 
ji^ for /i / 



v[p]^-i/{^p% F[p]^r// : (4.7) 



and we therefore introduce 



^o = -5p'-^, V[x{p)] = ^ + V[x{p)], (4.8) 
so that ( |4.5| ) is written as 

H^^i + V\x{p)]^i=Cj''il,i. (4.9) 



From Narnhofer ||2^ and following ||Tl|, we know that the operator Hq is not self-adjoint 
on the intersecting domain P(— (9^) n P(— l/(4p^)), with deficiency indices (1, 1). A family of 
self-adjoint extensions is determined by extending this domain to include the solutions of 

-dl[il^±^] - l/{Ap')^±, = ±iV±*, (4.10a) 

which in this case are 

= ^H^'^ (pe^-/4) , = V^i = VP^f (pe-*-/4) , (4.10b) 

where H^"^ and H^^ are the zeroth-order Hankel functions of the first and second kind. 

The real self-adjoint extensions are parameterized by the real angle r, and extend the oper- 
ator to act on the real functions 

:= i(e'"Vi + e-*"V'-i) = cos(T)3f?(^,) - sin(r)9(Vi), (4.11) 

where K and 9 represent the real and imaginary parts respectively. We write the extended 
operator Hq^t which acts like Hq on the extended domain 

ViHo,r) = V{-d^p) n P(-l/(4p2)) + {vi/^}, (4.12) 

and acts on '^r as 

Fo,,^, = - e-^"V-.) = - sin(T)5R(V'i) - cos(T)9(V'i). (4.13) 

We will require the following integrals: 



OO fCO 

2 / , /r^/ , N\2 



dp{^{'4)i)) = dp (9(V'»)) =: a « 0.1592, (4.14a) 
dp =: 6 « -0.1013, (4.14b) 



and using ( |4.1lD and ( [4.13| ) we find 

(^r^, ^^) = (a - 6sin(2T)) > 0, {^r, Ho^r^r) = -&cos(2t) G [-b, b]. (4.15) 

From the large-z behaviour of the Hankel function H^\z) ~ Y^^^6xp(i(z— 7r/4)) -|-o(z~'^/^) , 
we find the asymptotic behaviour of for large p 



'r ~ y^e"''/^ sin(/9/\/2 + vr/S + r) 



~ \l sin(/9/V2 + 7r/8 + r) , (4.16) 

and from the small-z behaviour of the Hankel function (7e is Euler's constant) Hq'\z) ^ 
1 + ^(ln(z/2) + 7e) + o{z'^ln{z)), we find near p=0 that 

*r ~ i cos(r) - ^ sin(r) - - sin(T) ln(p/2)) . (4.17) 
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At this point we note that only the r = extension (chosen imphcitly in |11|) corresponds to 
a perturbative field that is finite at p = and we will restrict ourselves to this case in what 
remains. 

From the form of the differential operator Hofi, for any functions tp,ip' £ I?(^o,o) we de- 
fine ip/3{p) := ^-nd it is straightforward to prove the scaling property: {^ip'p, Hofiipis) = 
P{ip' , Hofiip). We therefore consider the collection of 



(4.18) 



as the variational family, varying /3 to get an upper bound on the ground state energy. Using 
this scaling property and the results ( [4.14 ) and ( |4.15D , we find 



{^0,13, ^o,o^o,/3) = -bp, (I'd,/?, ^o,/3) = a/f3, 
and taking the expectation value of ( |4.9| ) leads to the bound 

c^' < /3{{^o,f3,V^o,p) - bp) /a =: ma^{0^). 



(4.19) 
(4.20) 



Naively we expect that very small values of /? will result in a bound of u)^ < /i^ since in these 
cases the bulk of the support of ^'o,/3 is shifted to large p and the asymptotic behaviour of the 



V (4/7) in the expectation value dominates the integral. In calculating the expectation value 
of the potential, we compute in practice 



dxe"2('^«-^o)y(x)*2 [p(^)]. 



(4.21) 



The effective potential is plotted in Figure ^ for some of the p = solutions, as is the depen- 
dence of max(u;^) on /3. We give the bound determined by minimizing max(a)^) with respect 
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Figure 4. The effective potential and behaviour of ma.x(u)'^) ( 4.20| ) with respect 
to P; the same values of s as in Figure |l] are given. 



to P in Figure 3(a). Although we will not be quantitative on this bound derived on the ^ 
solutions we have checked that it is consistent with the numerically determined u)^, generically 
providing a negative or zero bound. 
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5. Discussion 

In this work we have examined the Wyman sector of the nonsymmetric gravitational theory 
in some detail. We began by mapping it onto an Einstein-Klein Gordon system with metric 
g and scalar field tp with a positive semi-definite (coordinate-dependent) interaction potential. 
From this form we were able to employ a variational argument to show that there were neither 
globally regular solutions with nontrivial ip nor spacetimes with a horizon and nonzero ip in the 
exterior. Considering small r behaviour that matches the analytically known Wyman solution, 
we numerically solved the static background equations reproducing the analytically known 
Wyman solution and generating new solutions to the ^ field equations. 

The dynamical stability of these solutions was investigated by direct numerical determination 



of the lowest energy state and corresponding eigenvalue (similar to that appearing in |26|), as 



well as the variational approach of Jetzer and Scialom |11]. We found that the ^ = Wyman 
solutions are generically unstable and therefore not possible candidates for the endpoint of a 
collapsing NGT system, however they are most likely threshold solutions (no evidence of more 



than one unstable mode was found) and may therefore play a role in critical collapse [30|. 
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Appendix A. The Wyman Solution 

Here we develop some analytic results from the Wyman solution of the static, fi = 
field equations. Adopting the form of the solution given in jl^, we have (the coordinates are 
labelled by {t,a,6,(j)}; the functions 7, a, (3 and k are functions of the dimensionless radial 
coordinate a only) 

|g(Mi^)l = diag(7,-Q,-/?, -/3sin2(6')), gps] = k sin(6l), (A.la) 

where a G [0, 00], and we consider two separate solutions (corresponding to the ± signs through- 
out) 

7(a) =e±^ (A.lb) 
a{a) =4M|e^'"(l + s^)C{(t)-'^, (A.lc) 
f3{a) =8M|e^'^C((7)^^ (cosh(ap(7) cos(amO") — 1 -|- ssinh(ap(T) sin(am'7)), (A. Id) 

k((t) =8M|e^'^C((7)^^ (sinh(apO") sin(am(T) + s — scosh(ap(T) cos(am(T)). (A.le) 



Here 



(A.2) 



v 2 V 2 

s is a dimensionless constant of integration, and C{a) is defined by 

C(cr) := 2(cosh(apCr) — cos(amO")). (A--3) 

From its definition C(ct) and its derivative are non-negative, strictly increasing functions of u, 
and furthermore d(j[C{a)\ > C{a). 

We use the dimensionless coordinate x := r/(2Ms) where r is defined by requiring that the 
fundamental tensor be of the form ( |2.6|) , so that (3"^ + = and we find 

xl = e^^'/Cia), and d^[x] = -^x{d^[lnC{a)] ± l) < 0. (A.4) 
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For notational simplicity we have also introduced 

Xs:=x/{l + s'fl\ (A.5) 
Performing this coordinate transformation, we find the components of (|2.6|) 



[d„[C{a)]±C{a)) 

and if) may be determined from 

tan(V') - ^^'^^^"P^'' sin(amfT) + s(l - cosh(apfT) cos(a„(T)) 
ssinh(ap(T) sin(am(T) — (l — cosh(ap(7) cos(am(T)) 

For small a the relationship (|A.4| ) is easily inverted to give cr ~ l/x, which leads to C{a) ~ 
VT+~i2cr2 ~ and 

e'^^libl/x, e^~lTl/a^, ^"-5/(12x2), (A.8) 

justifying the identification of Mg as the asymptotic Schwarzschild mass parameter with the 
lower (upper) sign corresponding to the positive (negative) mass solution (in deriving the final 
result in ( [A.8| ) the numerator in (A. 7) must be expanded to fourth-order in a). These results 
are also consistent with the s ^ limit of (A.4), from which one finds = (1 it 1/x), and we 
find the negative and positive mass Schwarzschild solutions respectively (only the exterior in 
the latter case since the radial coordinate only extends to r = 2Ms in the limit). 
For large a we find Xg ~ exp(— (ap it l)a/2) , which may be inverted (s ^ 0) to give 

a ^ ln(x,), Cia) ~ e'^'"^ ~ ^-2a,/{a,±i)^ 

ap±l 

and from these we find 



[ap ± 1)^ 

Note that for the positive mass solution e'^ ^ as x ^ whereas for the negative mass 
solution e^ is singular as x — > 0; note that neither case satisfies the boundary conditions (|3.3| ) 
and are therefore not covered by the scaling argument of Section |3[ 

Deriving the small x behaviour of tp requires slightly more work. First note that da[ip] = 
—da[cos{ip)]/ sin(V') and we have 
2 

cos(^) = C~''^((7) (cosh(ap(T) cos(am'7) — 1 + s sinh(ap(T) sin(amO")) , (A. 11a) 

V 1 + 
2 

sin(-(/;) = C~^(o') (sinh(apcr) sin(ay„cr) + s — s cosh(apCr) cos(amCj)) , (A. lib) 

V 1 + 

from which we find in the large a limit that da[il^] ~ am- Integrating this and re-writing it in 
terms of x gives 

V ^ ln(x,) + tan-i(-s). (A.12) 

Op ± 1 

The constant is determined by noting that for a = an '■= mr/am, n = 1,2, . . . (n = Oisa 
special case since both the numerator and denominator of ( A. 7 ) become degenerate) we have 
from ( A. 7 ) that tan('0) = —s, and therefore 

'il^{an) =n-n + tan"^(-s). (A.13) 

As a check on the accuracy of the numerical solutions, the first six of these points are indicated 
by diamonds on the numerically generated solutions in Figure |l|. 
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Note that for this solution, the components of the metric g as defined in ( |2.8D are weh- 
behaved except at r = where g has a curvature singularity which shows up in the Ricci 
scalar (|2.9| ) as 

°^ -2a,/(a,-l)\ 

Furthermore, by considering radial null geodesies ending at r = we find that this is a strong 
curvature singularity. Radial null geodesies affinely parameterised by A are determined from 
the tangent vector {E is a constant of integration) u = [d\[t],dx[rf) = (^Ee~'^ ,±Ee~^'^~^^^^'^) 
which may be integrated near r = to find (choosing A = at r = 0, A increases with r and 
assuming that the trajectory passes through r = at t = 0) = (op±l)2Mst and r ~ 2Ms(H- 

sY^^{apEX/i2MsVTT^)Y^''^^'^^^^^^^- Using the result that R{u,u) = -E'^dr[e-^^+^^/r, 
we compute the derivative and write what remains in terms of the affine parameter to find 
X'^R{u, u) ~ a^/(2ap); the finite limit of this as A — > indicates that g has a strong curvature 
singularity at r = by an argument of Clarke and Krolak (see for instance []3l| ). 

Outgoing null geodesies are determined from dx[r] = ^e-('^+-^)/^ and upon integrating from 
r = (fixing A = at r = 0) and rewriting it in terms of a we find that 

pr poo 

/ = 2MsyrTs2 / daC{ay^ = EX. (A.15) 

Jo Jao 

This integral is strongly convergent as cq ^ oo, and has no poles except at o" = where the 
integral diverges. This indicates that the affine parameter A only goes to infinity as ctq ~^ 
(r — > oo); hence outgoing radial null geodesies reach any finite radius in a finite affine time, and 
the singularity is therefore visible (naked). Identical results may be derived from the so-called 
path equation (c./., |Q, Section 5, in particular (24)). 
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